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Adhesion between membranes is studied using a phenomenological model, where the inter- 
membrane distance is coupled to the concentration of sticker molecules on the membranes. The 
model applies to both adhesion of two flexible membranes and to adhesion of one flexible membrane 
onto a second membrane supported on a solid substrate. We mainly consider the case where the 
sticker molecules form bridges and adhere directly to both membranes. The calculated mean-field 
phase diagrams show an upward shift of the transition temperature indicating that the lateral phase 
separation in the membrane is enhanced due to the coupling effect. Hence the possibility of adhesion- 
induced lateral phase separation is predicted. For a particular choice of the parameters, the model 
exhibits a tricritical behavior. We also discuss the non-monotonous shape of the inter-membrane 
distance occurring when the lateral phase separation takes place. The inter-membrane distance 
relaxes to the bulk values with two symmetric overshoots. Adhesion mediated by other types of 
stickers is also considered. 
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I. INTRODUCTION 



Adhesion of membranes and vesicles has attracted considerable experimental and theoretical interest because of 
its prime importance to many bio-cellular processes . Theoretical treatments of membranes composed of single 
component lipid bilayers have revealed that generic interactions such as van der Waals, electrostatic or hydration 
interactions govern the adhesive properties of interacting membranes. It is also worthwhile to mention that related 
phenomena are found in unbinding transition of nearly flat membranes jsj or adhesion of vesicles to surfaces Q| . 

In addition to general non-specific interactions mentioned above, it is known from the works of Bell and coworkers 
|^,|| as well as others @], that highly specific molecular interactions play an essential role in biological adhesion. This 
interaction acts between complementary pairs of proteins such as ligand and receptor, or antibody and antigen. Well 
studied example of such coupled systems is the biotin-avidin complex. The avidin molecule has four biotin binding 
sites, two on each side, and forms a five-molecules biotin-avidin-biotin complex. The resulting specific interaction 
is highly local and short-ranged. Measurements by surface force apparatus or atomic force microscopy [p|jTc|] 
have shown that the force required to break a biotin-avidin bond is about 170pN. In related experiments measuring 
chemical equilibrium constants | pT[ , it was found that the biotin-avidin binding energy is about 30-35 ksT which is 
larger than thermal fluctuations. Other coupled systems are those of selectins and their sugar ligands where the bond 
is much weaker, of the order of bksT ^2 13|. 

More recently several models taking into account thermal fluctuations in membrane adhesion have been proposed. 
Zuckerman and Bruinsnia p^ , p^ used a statistical mechanics model which is mapped onto a two-dimensional Coulomb 
plasma with attractive interactions. They predicted an enhancement of the membrane adhesion due to thermal 
fluctuations. In another work, Lipowsky considered the adhesion of lipid membranes which includes anchored stickers, 
i.e., anchored molecules with adhesive segments |14 1^. It was shown that flexible membranes can adhere if the sticker 
concentration exceeds a certain threshold. If the multi-component membranes, including lipids and sticker molecules, 
undergo a phase separation, the adhesion is dominated by the sticker- rich domains. Further studies in this direction 
using mean-field theory and Monte Carlo simulations p^ ] obtained a phase separation which is driven both by 
attractive intra-membrane sticker interactions and fluctuation-induced interactions between stickers. 

The problem of multi-component membrane adhesion is intimately related to that of formation of domains (a 
lateral phase separation). This has been observed by several experiments. For example, the biotin-avidin interaction 
occurring during vesicle- vesicle adhesion was investigated by a micropipette technique ijl^ . The adhesion between 
one avidin-coated vesicle and a second biotinylated vesicle is followed by an accumulation of biotin-avidin complex 
in the contact zone. This accumulation of cross-bridges between the two vesicles is found to be a diffusion-controlled 
process. 

Adhesion-induced phase separation has been observed by Albersdorfer et al. and results from the interplay between 
long-range repulsive and short-range attractive interactions ||1^,q3. The membrane includes repeller molecules in 
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the form of lipopolymers (modified DOPE lipid with a polyethyleneoxide headgroup), mimicking glycocahx in real 
biological systems. The other component is a receptor molecule in the form of biotinylated lipids (DOPE-X-biotin). 
This lipopolymer is responsible for longer-range repulsive interaction, while the short-range attractive interaction is 
introduced by adding streptavidin to the extra-cellular solution. The streptavidin acts as a connector between the 
biotinylated lipids on the two membranes. A technique of reflection interference contrast microscopy pO| ] was used 
to observe domain formation on a vesicle adhering to a membrane supported on a solid substrate. The lateral phase 
separation on both membranes leads to the formation of domains of tight adhesion separated by domains of loose 
adhesion |l|,|9|. 

In a related work, adhesion between cationic vesicles and anionic supported membranes revealed that electrostatic 
interactions induce lateral charge segregation on the membrane [^,^ . This phase separation leads to patches of tight 
inter-membrane contact and decoupled "blisters" . Furthermore, adhesion of membranes including self-recognizing 
hemophilic molecules and lipopolymers has been investigated p3t . It was found that the initial weak adhesion is 
followed by slower aggregation into tightly bound domains coexisting with domains of weak adhesion. The result has 
been interpreted in terms of a double- well inter-membrane interaction potential due to the presence of the lipopolymers. 
Let us emphasize that in all the above mentioned experiments, it was reported that adhesion molecules aggregate 
spontaneously and form domains of tight adhesion. 

It is generally believed that multi-component biomembranes in physiological conditions are close to their critical 
point, and membrane functions are partially governed through phase separation processes. Moreover, concentration 
fluctuations in the vicinity of the critical point may affect biophysical properties of membranes and can be of impor- 
tance in regulating membrane processes in a robust way. Recently this conjecture was supported by an experiment 
of an insoluble Langmuir monolayer at the air/ water interface ||2^,^^. The monolayer was prepared in two different 
steps. The first mimics the composition of the inner leaflet of a cell biomembrane, while the second mimics the outer 
leaflet. In both cases, by using fluorescence microscopy technique, it was found that the Langmuir monolayer is close 
to its corresponding critical point of demixing. 

So far, the interplay between lateral phase separation and membrane adhesion has not been considered theoretically 
in detail except in Refs. |l4|-|l^]. The work in Refs. pl| , p2| deal only with the specific case of oppositely charged 
membranes. In this paper we provide a general phcnomcnological approach for the adhesion of multicomponent 
membranes. Using a mean-field theory, we investigate how the lateral phase separation within the membrane is 
affected by the adhesion of membranes. Like in Refs. p^-p^, we consider adhesion mediated by sticker molecules. 
Sticker molecules are polymers or macromolecules anchored to one membrane and interacting with the other membrane 
by another sticky part of the molecule. They can form bridges between two adjacent membranes (so called trans- 
interaction) 10, and play an essential role in the adhesion of cell membranes in biological systems. 

We distinguish three types of adhesion depending on the structure of bridges as represented in Fig. P. (i) "Bolaform- 
sticker" adhesion where each bridge molecule consists of a single sticker having two sticky ends (Fig. |l|(a)). One sticker 
end is anchored to one membrane while the other end is adhering directly to the second membrane, (ii) "Homophilic- 
sticker" adhesion where the bridges are formed by two stickers of the same type (Fig. |l|(b)). Each sticker is anchored 
on one of the membranes, while their free ends bind together to form the bridge, (iii) "Lock-and-key" adhesion where 
the bridges consist of two different stickers forming a ligand-receptor type bond (Fig. |l](c)). This case represents 
an asymmetric adhesion due to the lack of symmetry between the ligand and receptor. In the present work, we 
mainly discuss the symmetric bolaform-sticker adhesion (case (i) above) using a model where the equilibrium spacing 
between two membranes is coupled to the local concentration of stickers. Even in the latter symmetric case, a certain 
asymmetry can be obtained by controlling separately the sticker chemical potentials on the two membranes. An 
important consequence of our model is that the lateral phase separation is enhanced. 

This paper is organized as follows. In the next section, we explain our phenomenological model of bolaform-sticker 
adhesion. The mean-field phase diagrams are given in Sec. III. The inter- membrane distance between two coexisting 
domains is calculated in Sec. IV. Finally discussion is provided in Sec. ^ where the other types of adhesion mentioned 
above are considered. 



II. BOLAFORM-STICKER ADHESION 



In this section, we treat the case where the adhesion is mediated by a single type of sticker molecules which are 
anchored irreversibly to one membrane and stick to the other membrane by another sticky part of the molecules as 
in Fig. |l|(a). The anchor segments consist of a hydrophobic segment and penetrate into the hydrophobic interior 
of the lipid bilayer. The sticky segments, on the other hand, adhere directly to another membrane having some 
potential of sticking As mentioned above, we call this a "bolaform-sticker". Consider two interacting membranes 
labeled by i = 1,2 consisting of lipid molecules and bolaform-stickers as schematically shown in Fig. 0. Let the 
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sticker concentration in each membrane be denoted by ^i{r), where r = [x, y) is a two-dimensional planar vector and 
< i^iir) < 1. Note that the average concentrations of stickers on the two membranes, {i/^i) and (^'2)7 do not have 
to be the same. 

When the adhesion molecules are very flexible, they can bend back to form arches on a single membrane. In order 
to avoid such a situation in experiments and in the model, the bending rigidity of the sticker molecules should be 
sufficiently large. Hence we assume that for stiff enough stickers all the bonds are intcr-mcmbrane ones connecting 
the two separate membranes as considered in Refs. JT^ , p^ . 

The interaction between two stickers on the same membrane is called cis-interaction and can be repulsive or 
attractive. Here we discuss the case in which this interaction is attractive. Then, below a certain critical temperature, 
the multi-component membrane undergoes a first-order phase transition and stickers form lateral domains. As shown 
in Fig. ^ a sticker-poor phase coexists with a sticker-rich phase in the two-phase region of the phase diagram. The 
sticker critical concentration ipc and the critical temperature Tc are assumed to be the same for the two planar 
membranes. We define the concentration difference (/),;(r) for each of the membranes with respect to the critical 
concentration by 



(2.1) 



where varies between ^1 < 4>i < 1- 

The total free energy of the two coupled membranes can be written as a sum of several terms detailed below. 
The first contribution describes the lateral phase separation of each membrane. Motivated by recent experiments 
on Langmuir monolayers |^,^ demonstrating that the inner and outer leaflets of biomembranes are close to their 
critical point, we employ a phenomenological Ginzburg-Landau free energy which is an expansion in powers of the 
order parameters Hence we have 
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(2.2) 



This expansion for the free energy can be justified close to a critical point where the (/i.^'s are small enough. The 
parameter c representing the line tension acting at the domain boundary, and the reduced temperature t — {T~Tc)/Tc 
are taken to be the same for the two membranes. On the other hand, the chemical potential fii, coupled to the 
membrane sticker concentration 0i, can differ between the two membranes since the sticker concentrations on the 
two membranes do not have to be the same. We recall that each bolaform-sticker is modeled with one of its ends 
anchored irreversibly to one membrane, while the second sticky end is attracted to the second mem bran e. The sticker 
concentration is associated with the anchored end of the stickers. The normalization factor 1/2 in (2^) is introduced 
in order to write down the free energy per single membrane. The coefficient of the fourth order term can generally 
be set as a positive constant without loss of generality. 

It is convenient to introduce the following new variables for the average and the difference between the two concen- 
trations: 



62 +01 ■02 + "01 



V'c, 0- 



01 02 - V'l 



(2.3) 



where — 1 < 0+ < 1 and — 1 < 0_ < 1. In terms of these new variables, (2.2) can be written as 



dV 



ic[(V0+)2-H(V0_)2]-Hii(0^-H0^_) 



+-{(j)X + 60^02_ + 0i) - M+0+ - M-0- 



(2.4) 



where 



Ai2 + Ail A*2 
7^ ' Ai- = 



Ml 



(2.5) 



The chemical potential /i_ associated with the order parameter 0_ is non-zero when the symmetry between the two 
membranes is explicitly broken. Namely, the two interacting membranes have different average concentrations of 
stickers. 

Next we consider the out-of-plane deformation energy of the two membranes. As depicted in Fig. ^, the membrane 
shape is parameterized by their heights ii{r), t2{r'), above the x-y reference plane. Working in the Monge represen- 
tation it is implicitly assumed that the membranes remain flat on average and have no overhangs. This approach can 
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be also useful to treat adhesion of vesicles in their contact zone. When the vesicle is large enough, it will be roughly 
flat close to the contact region, and the entire vesicle can be thought of as a reservoir for the stickers. Returning to 
the deformation energy, it can be written as the sum of the bending energy and the surface tension of each of the two 
membranes separately, as well as the interacting potential energy between them |26|,|27|: 



dV 



+ / dVi-(£i-£2;V'i,V'2), 



(2.6) 



where n is the bending rigidity, a is the mechanical surface tension acting on the membranes, and v is the potential 
energy per unit area representing the intcr-membrane interactions. For simplicity, k and a are assumed to be equal for 
the two membranes and do not vary as a function of the sticker concentration ipi. The potential v{ii — ^2; V'l: V'2) can 
be generally assumed to be a function of the local relative height coordinate £1 — £2 and the sticker concentration ipi. 
The former assumption is the so-called Derjaguin approximation The dependence on the sticker concentration 
tpi will be considered later. 

We now make a change of variables and transform to the center-of-mass and relative coordinates given, respectively, 
by 



L = 



(2.7) 



Only terms which depend on £ can be considered in the case where the center-of-mass is stationary, hence L is a 
constant of motion. Then (2.6) can be written as 



(2.8) 



In the high-tcmperature phase, the stickers are homogeneously distributed, and each of the membrane is in a one- 
phase region on the phase diagram. We assume that even in the absence of sticker molecules, such membranes are 
bound to each other due to the balance between the short-range repulsive {e.g., hydration interaction) and longer- 
range attractive interactions {e.g., van der Waals interaction). Hence, we do not consider the interesting problem of 
the unbinding transition |14ilq|. Although the membranes are always bound together, their equilibrium distance 



depends on the sticker concentration. Let us consider the potential v{t, -iAii "02) in (^-^ ) for 0+ — 0. Note that 0+ = 



means that ("01 + V-'2)/2 — ?Ac, namely, the average sticker concentration on the two membranes is at its critical value 
^c- The inter-membrane potential v{£] 0+ — 0) is assumed to have a single minimum at a certain inter-membrane 



distance 



€0 for 



0. This gives the equilibrium distance between the two bound membranes. The deviation of 



the inter-membrane distance from is defined by the dimcnsionless quantity 5{r) given by 



5{r) 



£{r) 
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(2.9) 



For small deviations from the minimum of the potential, v{£; 0+ = 0) can be expanded to second order. This is known 
as the harmonic approximation and gives 



v{£;cf>+ = 0)^v{£o) + ^v"{£o){£-lof 



(2.10) 



where u"(€o) is the se cond de riva tive of v with respect to £ evaluated a,t £ = £0, V = v"{£o)£q, and v'(£) = at = £0 
Using (jU) and ( ^0|) , dJ) can be written as 



F2 



dV 



(2.11) 



with K = k£q and S 



This is the expression of the deformation energy within the harmonic approximation and 
it served as a starting point to many calculations on membrane adhesion . 

Now we will include the effect of the adhesion on the phase separation and suggest a lowest-order coupling between 
the composition (j3i(r) and the inter-membrane distance d{r). When the membranes are quenched into a two-phase 
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region of the phase diagram a sticker-poor phase coexists with a sticker-rich phase. As shown in Fig. this can lead 
to different inter-membrane distance for the different membrane domains. Since the sticky segments of the bridges 
adhere directly onto the two membranes, the coupling is proportional to the sum of the local sticker concentrations 
of the two membranes. This can be phcnomcnologically represented by the following coupling term: 



(2.12) 



where the coupling constant a is positive preferring smaller sep arati on (5 < in regions where the average concentration 
(f)+ is positive (or i)! + "4)2 > 2ipc)- In the last expression of (2.12), we have neglected the linear terms in and 6, 
which merely shift the chemical potential or minimum of the potential, respectively. Depending on the value of 
this coupling term not only introduces a shift of the minimum of the potential v but also changes the minimum value 
of the potential. Note also that ( 2.12 ) is symmetric with respect to the exchange of the two membranes 1^2. 

The above linear coupling energy can also be understood in the following way. Let us first consider a single flexible 
membrane with sticker molecules adhering to a flat substrate. Suppose Vr{i) and Vf{i ) are the potentials for sticker- 
rich and sticker-free membrane, respectively. Following the same discussion as in ( [2.10 ), each of the potential is 
parabolic around a different separation: Vr « a{£ — l^Y and w/ ~ a{i — (-f Y ■ The effective potential can be obtained 
by a linear combination of these two potentials, i.e., (1 — ip)vf{i) + '0Wr(^)- By expanding Vf{£) and Vr{£), we get a 
coupling term which is proportional to ipf. In the case of adhesion betwee n two membranes, we add the contributions 
from both of the membranes and obtain the coupling energy as given in ( 2.12| ). The same argument can be repeated 
for any arbitrary adhesion potentials, Vf{£) and Vr{£), provided each of them has a single well-defined minimum at 
some distance £. 

The total free energy considered in our model is the sum of (|2j), ( pHI ), and ( ^l2| ): 



F = Fi + F2 + F3. 

Here it is convenient to convert to Fourier space. The Fourier transform of any function /(r) is defined as 

f{q)= I d\f{r)e^1-^, 



(2.13) 



(2.14) 



where q is the two-dimensional in-plane wavevector. The total free energy can be expressed as 

'1 



F = 



-y 



;(t + cg2)(|<^+(9)|2 + |^_(g)r 



dV 



-(0 



(2.15) 



where A is the area of the membranes projected on the x-y plane. For convenience the free energy ( 2.15|) is written 
as a combination of real space and Fourier space terms. 

Within the mean-field level, the free energy in terms of 0+ and is obtained by functionally minimizing F with 
respect to <5(q). Then we find 



a4>+{q) 



V + J:q^+ Kq- 



A ■ 



(2.16) 



Hence the inter-membrane distance (5 = {£ — £13) I £q is fully determined by the value of 0+. By inserting ( 2.16 ) into 
(2.15), the resulting free energy depends only on and and becomes 



F = 



if+(9)|0+(9)P + i(i + cg2)|^_(g)|2 



dV 



+ 60^0^ +<^i) 



(2.17) 
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where 



If we expand the last term in ( ^.18 ) for small q, we obtain 



T+iq)^it-j)+[c+^)q\ (2.19) 



with 



7^—. (2.20) 



The parameter 7 is an important parameter characterizing the coupling strength. The first two terms in ( 2.19| ) implies 
an upward shift of the transition temperature, as will be discussed in detail in the next section. We also find that the 
presence of the coupling {a 7^ 0) increases the line tension c provided the mechanical surface tension S is non-zero; 

III. PHASE DIAGRAMS 

In this Section, we calculate the mean-field phase diagrams for bolaform-sticker adhesion using the free energy 
explained in the previous section. In order to study the b ulk pr operties of the system, we set q = and study the 
homogeneous solutions, (pi^s and S being constants. From ( |2.16|) , the inter-membrane distance which minimizes the 
free energy is given by 

(3.1) 



V 



Since a is positive, S is negative (smaller inter-membrane distance) for positive 0+, and S is positive (larger inter- 
membrane distance) for negative . By substituting back this S into the free energy / per unit area for homogeneous 
(constant) and <p-, we obtain 



f = lit -1)^1 + ^tc^'- 
1 



^ +6<j>icj)i+cbt) -(!+<!>+ -fi-cj)., (3.2) 



where 7 is defined in ( ^.2C ). Notice that 7 is never negative and vanishes only when a — 0. Therefore, when 



/i+ = /i_ = 0, the field will order before and the phase with 7^ and 0_ = is expected Although 
the phase behavior of this free energy can be examined in general, we concentrate here only on two particular cuts in 
the parameter space, i.e., ji- = and ^+ = 0. In these cases one can clearly see the effect of adhesion on the lateral 
phase separation. 

A. The case /i_ = 

When /i_ — the two membranes have the same chemical potential /ii = ^2- Since the chemical potential is 
coupled to 4>+ , / can be minimized first with respect to . A "symmetric phase" is obtained for t + 30^ > with 

= 0, (3.3) 

where the two membranes have the same concentrations, (j)! — 4'2- Likewise, two "asymmetric phases" are obtained 
for t + 30+ < with 



=±.^-<-30^. (3.4) 
In the asymmetric phase the two membranes have different concentrations <j)i ^ <i>2. After inserting these expressions 



into (3.2) with ^_ = 0, the free energy becomes 
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/l 



for 
for 



t + 30^ > 



t + 301 < 



(3.5) 



Notice that this free energy is continuous at i = —30'^. This free energy fi can now be minimized with respect to 
. The resulting equation of state is written as 



AH = 



{t - 7) 0+ + 
- (2i + 7) 0^ 



for i + 30| > 
for i + 30| < 



(3.6) 



The phase diagram can now be calculated and the two-phase region is obtained by the Maxwell construction. The 
phase diagram for = is illustrated in Fig. ||. 

For t > 0, only the symmetric phase with 0_ = can appear since t + 3<^^ > 0. Two symmetric phases with 
different can coexist when t < "f. The coexistence curve is simply given by 



and the associated critical point is located at 



(<,(/.+ ,Ai+)c = (7,0,0). 



(3.7) 



(3.8) 



We stress that the critical temperature is increased from tc — to tc = 7 = a^/V due to the coupling between 
the composition 0+ and the inter- membrane distance S as given in (2.12). In other words, the phase separation is 
enhanced by the adhesion of membranes. As prese nted in Fig. ^, the two coexisting values of 0+ given by (3.7) lead to 



different inter-membrane distances 6 according to (3.1). Since a > 0, (5 is negative {£ < £0) in the sticker-rich domain, 
and this phase is called the "tight phase" (T). On the other hand, 5 is positive {£ > £0) in the sticker-poor domain 
and this phase is called the "loose phase" (L). However, for each of the coexisting tight and loose phases, 0_ = 0, 
which means that the sticker concentration is the same in the two membranes, 0i =02- 

For —7/2 < t < 0, th e as ymmetric phase with ^ is always unstable, and the tight and loose phases coexist 
according to (3T) and (3/7). For t < —7/2, the asymmetric phase can be locally stable but it is only metastable. 
Namely, its free e ner gy is higher than that of the symmetric phase. Hence the coexistence between the tight and loose 
phases given by (3.1) and (3.7) preempts the asymmetric phase. The limit of metastability of the asymmetric phase 
is obtained by calculating the second derivative of the second equation of ([3.5D with respect to (p^. This leads to 



± 



-2i-7 



24 



(3.9) 



which is also shown as a dotted line inside the L-f-T coexisting region of Fig. ||(a). 

In summary, for = 0, the asymmetric phase 7^ does not exist as a stable phase for any temperature. 
At most it is metastable and occurs within the L-l-T coexistence region. The tight and loose phases coexist for 



t < 7 = a'^/V according to (3.1) and (3.7) 



B. The case ^+ = 



Next we consider the case of /i+ = but with fj,^ 7^ 0. This means that the chemical potentials of the two 
membranes have the same magnitude but opposite sign, i.e., fii — —^2- This is a special c ase of the more general 
situation where the symmetry between the two membranes is explicitly broken. Now / in (|3.2| ) can be minimized 
with respect to (f>j^ first. As long as t-\- i(fpL > 7, the only solution is 



=0. 

This is called the "middle phase" (M) where the inter- membrane distance is exactly £0 (or S — 0). 
(/>+ = means that ipi + 1IJ2 = 2^/'c- For t + 30?_ < 7, we have the tight (or loose) phase with 



(3.10) 

Again note that 



± 



-i + 7 - 3(/)2 



where £ deviates from £0 (or 6^0) according to (3.1). Since /i+ 
energetically degenerated and they coexist. By substituting back into (3.2) with /i 



(3.11) 

0, both the tight and the loose phases are 
0, the free energy becomes 
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^ / + \(t>t - li-<t>- for t + 3<^2_ > ^ 

1 -i(t-7)' + i(-2t + 37),/)^-20i-//_0_ for t + 3</)2_ < 7 ' ^''■^^^ 

After minimizing with respect to (p- , the equation of state is given as 

_ / i0_ + for t + 30!, > 7 , . 

" \ (-2i + 37) 0_ - 8(/.^ for t + 30^<7- ^'^■^'^^ 

The calculated phase diagrams for — are shown in Fig. ^. The phase diagram is symmetric about 0_ = 
and = as a consequence of the coupling term, and lack of any odd terms in 0+ in the free energy. For 

t > 7 there is a one-phase region of the middle phase with 0+ = since t + 30^ > 7. For 57/6 < f < 7, the system 
undergoes a second-order phase transition between the middle phase (0+ = 0) and the tight (or loose) phase (0+ ^ 0). 
The analytical expressions of the second-order phase transition lines are 



±a/-^, (3.14) 



and 



2i + 7 /^i + 7 



in Fig. ^, respectively. On the second-order phase transition line, 0+ goes continuously to zero. 

For t < 57/6, the transition changes to first order. This has been numerically determined by the Maxwell construc- 
tion. The point which connects the first- and second-order phase transition lines is a tricritical point p^ ]. In our 
model, it is located at 

itA-,^^-)te, = (^7, ± ^7^/^ ± ^7^/^) . (3.16) 

The first-order phase transition corresponds to the coexistence of the middle phase with 0+ = and the tight (or 
loose) phase with 0+ ^ 0. The obtained two-phase coexistence region is indicated by "M-l-T" in Fig. ^(a). Within 
the present Ginzburg-Landau expansion, the tight phase persists even if we go to low temperatures. Because of 
the degeneracy between tight and loose phases, the first-order line near the tricritical point actually corresponds to 
coexistence of three phases: tight, loose, and middle phases. 

In continuation to the discussion of the previous subsection (/i_ = 0), we see that the phase separation is also 
enhanced for the = parameter space. It occurs at higher temperatures, since the tricritical temperature Ucp = 
57/6 — 5a^/6V is positive for a ^ 0. It is important to notice that in the middle phase with 0+ — 0, the inter- 
membrane distance is £q since (5 = 0. On the other hand, in the tight (loose) phase with 0+ > (0+ < 0), according 



to (3.1) and ( |3.7[ ), £ < £0 {i > £o)- The coexisting membrane domains between tight and middle phases, or between 
loose and middle phases is schematically represented in Fig. 0. 

We end this section by commenting on the general case when both /i+ and /i_ are non-zero. When becomes 
non-zero, the degeneracy between the tight and the loose phases is lifted. In such a case, instead of the three-phase 
coexistence for /i+ — 0, there is a coexistence between either tight and middle phases, or between loose and middle 
phases as shown in Fig. ^(a) and (b), respectively. Notice that the tricritical point exists only when fj,+ — 0. In a 
more general phase diagram drawn in the {t, 11+, space, three second-order lines meet at the tricritical point. In 
the three-dimensional parameter space, these second-order lines lie on the perimeter of two-phase coexistence planes 
between either tight and middle phases (T-l-M) or between loose and middle phases (L-l-M). 



IV. NON-MONOTONOUS MEMBRANE PROFILE 



One of our assumptions was that the inter-membranes potential v{£] = 0) has a single minimum aX i ~ £q when 
0_i_ = 0. In the absence of thermal fluctuations, two homogeneous membranes are bound with inter-membrane distance 
£q for 4>+ = 0. In this section, we calculate the profile of the inter-membrane distance between two membranes which 
are quenched below the phase separation temperature. 

We first expand the potential v{£] 0+ = 0) up to the fourth order terms in (5; 

v{£-<l,+ =Q)^v{£^) + ]^V6^ + hj5\ (4.1) 
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where V = w"(^o)^o before and U = v^^\£o)£q > 0. Both V and U axe positive constants because of the convexity 
of V at its minimum. Suppose that each of the mem bran es is in its high-temperature phase (t > 0) . Then the fourth- 
order (f)i terms in the Ginzburg-Landau expansion (2.4) can be neglected since the second-order terms are positive. 
The resulting free energy with /i+ = = is 



F 



dV 



1 



4 



-y 



i(i + cg2)(|^+(g)p-h|<^-('7)P) 



dV -ub\ 



We now minimize F with respect to the concentrations (f>+{q) and (j)^{q) and obtain 



t + cq- 



:~S{q), 0-(q) = O. 



By inserting these equations into (4.2) and expanding for small q, the free energy can be written as 

1. 



F ■ 



dV 



iife(V2<5)2 -I- lEe(V<5)2 + ^VJ^ + -U,6' 



(4.2) 



(4.3) 



(4.4) 



with 



K, = K- 
V, = V- 



t 



u, = u. 



t2 



(4.5) 



We see that for i > the coupling always increases the mechanical tension So > S, but reduces the rigidity < K 
and the potential strength 14 < V. 

Let us consider the strong coupling case when Vc <0 but still having > 0, namely, 



ft 

V < — < K { - 

t \ c 



(4.6) 



For t > 0, although no phase separation occurs in the absence of coupling (a = 0), it occurs for non-zero a. The 
minimum free energy of the membranes is given by solving the Euler-Lagrange equation obtained by minimizing (4.4) 
with respect to the inter-membrane distance S: 



The two uniform (bulk) solutions of (4.7) are 



(4.7) 



(4.8) 



We assume a one-dimensional profile 5{x) describing the inter- membrane distance along the x-direction. A typical 
profile determined by a numerical solution of the Euler-Lagrange equation (4.7) using a relaxational method is shown 

in Fig. 1^. It is convenient to rescale the variables 5 and a; as C = Uc^^S and u = Kc ^^'^ul^^'^x, respectively, yielding 
the following one-dimensional profile equation: 



c + c 



(4.9) 
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Only two independent combinations of the four parameters K^. Se, T4 and Uc exist. In Fig. ^ they are set to be 
Yiel {kI^'^uI^^) = 0.1 and Vc/Uc^^ — —1, respectively. The profile has a large slope at the interface x = u = 0, 
but relaxes to the bulk values ±(5o at a; = ±oo in a non-monotonic fashion with two symmetric overshoots, having a 
height greater than Sq. These overshoots are suppressed by increasing Se or by increasing the coupling strength a. 
The maximum value of S at the overshoot scales as as can be seen from ( |4.8| ). The overshoot of the profile 

is followed by a damped oscillation which minimizes the curvature energy. This behavior is similar to the nonlinear 
response of membranes to local pinning sites ||3^|35[| or membranes adhering to a geometrically structured substrate 
p9| and is a result of the 4th order derivative in the profile equation. The oscillatory decay has been also predicted 
for the membrane profile between two inclusions such as proteins p6|^-p9|. 

The configuration of the phase separated membranes corresponding to the above inter-membrane distance S is 
schematically represented in Fig. ||. In the case of the adhesion of a single flexible membrane onto a supported 
membrane, the supported membrane cannot have any shape fluctuations. Therefore, the inter-membrane distance 
profile calculated in this section can be regarded as a distance of the flexible membrane from the substrate with 
respect to its equilibrium distance £o- 



V. DISCUSSION 



A. Main Findings 

In this paper, the interplay between adhesion and lateral phase separation of multicomponent membranes is investi- 
gated. We consider the "bolaform-sticker" adhesion where adhesive bridges are formed by a single sticker having two 
sticky segments and adhere directly onto the two membranes, as shown in Fig. |^(a). We proposed a phenomenological 
free energy cons istin g of three parts: (i) the free energy describing the lateral phase separation of stickers on each 
membrane (see (^^)); (ii) the deformation energy of the two membranes, which is the sum of the bending energy, 
the surface tension, and the potential energy (see (|2.6|)); and, (iii) the coupli ng en ergy between the inter-membrane 



distance and the average concentration of stickers on both membranes (see ( 2.12 )). The difference of the chemical 
potentials between the two membranes is also taken into account because the sticker concentrations do not have to 
be the same. 

We calculate the phase diagrams describing the bulk properties for two particular choices of the chemical potentials, 
i.e., fi- = (/xi = ^2) and /i+ = (/xi = — /i2)- In the c ase o f /i_ = 0, the critical temperature increases depending 



on the coupling strength and the potential strength (see (3_^)). Hence the lateral phase separation is enhanced due 



to the adhesion. This is one of the main consequences of our model. When the phase separation takes place, the 
inter- membrane distance is smaller for the domains rich in the sticker molecules ("tight phase"), and larger for the 
domains poor in the stickers ("loose phase"). In the case of fi+ = 0, our model exhibits a tricritical behavior. The 
upward shift of the tricritical temperature also indicates the enhancement of the lateral phase separation. 

We find that the line tension for the lateral phase separation increases because of the coupling effect as long as the 
mechanical surface tension is non-zero. We have also calculated the inter-membrane distance profile between the two 
membranes which are quenched below their phase separation temperature. Because the membrane shape is governed 
by the bending rigidity, the inter-membrane distance profile relaxes to the bulk values in a non-monotonic way with 
two symmetric overshoots. 



B. Membrane Adhesion on Solid Surfaces and Supported Membranes 

So far, we have mainly discussed the adhesion of two membranes. Our model also applies to the case where a 
single flexible membrane with sticker molecules adheres to a flat substrate or a supported membrane Let us 

discuss these two cases separately. For a flat substrate without any supported membrane on it , the contributions 
from the second membrane (say i = 2) can be dropped from the model. The coupling term ( ^.12| ) simply reduces to 
a0i(5 because the stickers are assumed to adhere directly to the substrate. The second case is that of a supported 
membrane with sticker molecules. Unlike the case of two fluctuating membranes discussed in Sec. |8[ the supported 
membrane does not have any shape fluctuations. However, even in such situations, there is an enhan cem ent of the 
phase separation due to the coupling effect and the upward shift of the critical temperature is given by ( [3.8[ ) . Another 
related situation is the case where a membrane is composed of two different lipids and the membrane is put close to 
a flat substrate. If the two lipids feel different hydration force and prefer different distances from the substrate, the 
phase separation between the two components will be enhanced by the adhesion for the same reason described in this 
paper. 
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C. Relation to Other Models 



There exists an analogy between the phase behavior of our membrane system with that of metamagnets (magnets 
which undergo fist-order phase transitions in an increasing magnetic field) or ^He-^He mixtures described by the BEG 
(Blume- Emery-Griffith) spin-one model |^,^. Moreover, the phase diagrams for = are analogous to those 
describing the phase separation of two-component mixtures in fluid bilayers which also exhibits tricritical behavior 
p2| and other related amphiphilic systems ]4^-|4^. However, in the former case of two-component bilayers the 
concentration difference between the two leaflets of the membrane is linearly coupled to the curvature of the bilayer 
and the difference in the chemical potential is not taken into account. 

In our paper, we did not address the problem of the unbinding transition. We rather assumed that the membranes 
are always bound together, even in the absence of any sticker molecules. This assumption is partially motivated by 
the experimental study of Ref. |2^] where suspended membrane (part of the giant vesicle) was claimed to be bound 
to the supporting membrane even in the absence of sticker molecules. In this case, the inter-membrane distance t 
stays finite and it is permissible to expand the free energy around the minimum. Hence the phase separation consists 
of loosely and tightly bound patches. The interplay between unbinding transition and phase separation of multi- 
component membranes has been considered in other theoretical works p^-p^. The adhesion there is only brought 
about by sticker molecules, and the phase separation is induced both by attractive interactions and fluctuation- 
induced interactions between the stickers. Although their model treats a different aspect of the more general problem, 



the fluctuation effect yields similar consequences compared to ours. We assumed that the cis-interaction in (2.2) is 
attractive, and tracing over the inter-membrane distance 5 yields a term proportional to ^(\)\ ■ Since this term does not 
depend on the sign of 0+ , it has a similar effect as fluctuations although our treatment is restricted to the mean-field 
level. 

It is worthwhile to comment here the difference between the present study and that of Ref. j2^. In their paper, it 
is found that the adhesion between the membranes including homophilic recognition molecules and repeller molecules 
is controlled by lateral phase separation. The multiple competing states of adhesion is attributed to the double- 
well inter-membrane interaction potential generated by the competition of two forces; attraction between homophilic 
molecules and the repulsion between repeller molecules. By changing the repeller concentration, the double-minimum 
potential causes the first-order transition between a state with inter-membrane spacing set by the thickness of the 
repeller molecules to a state with a spacing set by the bare potential (van der Waals plus hydration interactions). 
In our work, the effect of repeller molecules is not ta ken in to account and the minimum of the potential depends on 



the sticker concentration through the coupling term ( 2.12 ). When the stickers are phase separated and two different 



values of the sticker concentration coexist, the inter-membrane potential has double-minimum. However the physical 
origin of this double-minimum potential is different from that in Ref. because it is not due to the presence of 
repeller molecules. 

D. Other Types of Sticker Molecules 

As mentioned in the introduction, "honiophilic-sticker" adhesion occurs when the adhesive bridges are formed by 
two stickers of the same type bound together by their two sticky segments (see Fig. 0(b)). Suppose tpi (z = 1,2) 
denotes the sticker concentration on each membrane. Then the inter-membrane distance depends on the product of 
each sticker concentration expressing the probability to have two stickers - one on each membrane - at the same 



position. Using (2.1) this coupling term can be written in terms of 0^ as 

ll^li>2^<f>l(t)2+1pci(t)l+(t)2)+1pl- (5.1) 

An interesting remark can be made for homophilic-sticker adhesion. The resulting phase separation within each 
membrane leads to three different values for the inter-membrane distance. The inter-membrane distance between 
domains rich in stickers on both membranes (rich-rich), as well as between rich-poor domains, and poor-poor domains 
can be different |^,^. Notice that these three different inter-membrane distances correspond to the tight, middle, 
and loose phases in our model. 

A third case is that of "lock-and-key" adhesion due to the formation of chemical bonds between lock-and-key 
types of stickers, e.g., ligands and receptors (see Fig. 0(c)). Suppose that both types of stickers are distributed on 
the two membranes and V'i now represents the local concentration, say, of the lock molecules. First let us assume 
that the membranes are saturated with sticker molecules (no lipid). Then, 1 — represents the concentration 

of key molecules. Since domains rich in lock (key) molecules on one membrane adhere with domains rich in key 
(lock) molecules on the other membrane, the coupling tern in the free energy Ft, will be a coupling between the 
inter-membrane distance and 
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V'l(l-V'2) + (l-^l)V'2 

-201,^2 + (1 - 2^e)(0i + 02) + 2V;,(1 - Vc). (5.2) 



This term is symmetric with respect to the exchange of two membranes. Both in (xl) and (|5.2D , we see that the 
lowest order term in the concentration (except the constant term) is proportional to (/)+. If there is a linear coupling 
between the inter-membrane distance and (/>+ in these cases, we expect an upward shift in the transition temperature 



and the phase separation will be enhanced as argued above. Due to the presence of higher order terms in (p.lD and 



(5.2), however, the phase behavior will be more complex. 

Let us now take into account the presence of lipids in the lock-and-key adhesion. If a single type of sticker is present 
on each membrane, namely, lock molecules on membrane 1 and key molecules on membrane 2, we can regard ■01 
and ^02 as the concentrations of lock and key molecules embedded in the lipid membrane, respectively. Then 1 — ^/^i 
and 1 — "02 describe the concentration of the second component (lipid) on each membrane, respectively. In this case. 



the inter-membrane distance depends on ^i'02 as in (5.1). When both lock and key molecules are present on both 



membranes, one has to start with a three-component mixture for each of the membranes. The generic lattice model 
to study the behavior of ternary membranes of monolayers is the BEG spin-one model |^,^,^ . Here one has to 
include the coupling between the two membranes. If we denote the concentration of lock and key molecules on each 
membrane as 0^ and t/;/^ (and hence the concentration of the dilu tion lipid is 1 — T/^f" — V"/^); the inter- membrane 
distance now depends on '0f 02^ + '01^02' which is similar to (5.2). More detailed calculations for the homophilic 



stickers and lock-and-key stickers and their influence on membrane adhesion are left for future studies. 
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FIG. 1. Schematic representation of various types of adhesion between two membranes. The two membranes are represented 
by black hnes. (a) Bolaform-sticker adhesion: bridges consist of a single type of sticker molecules which are anchored to 
one membrane (filled circle) and stick to the other membrane by another sticky part of the molecules (open circle), (b) 
Homophilic-sticker adhesion: bridges consist of two identical stickers which are bound together by their respective sticky end 
segments, (c) Lock-and-key adhesion: bridges consist of two different stickers forming a ligand-receptor type bond. 

FIG. 2. Two adhering membranes. The reference x-y plane is shown as a dashed line. The height of two membranes 
measured from this plane is denoted as ii and £2, respectively. The sticker concentration on each membrane is denoted by tpi 
and ip2, respectively. 

FIG. 3. Schematic phase diagram for a single membrane containing sticker molecules. The concentration of the sticker 
molecule is tp. The critical concentration and the critical temperature is denoted by and Tc, respectively. Within the 
coexistence curve the membrane separates into sticker-rich and sticker-poor region (A-l-B coexistence). The membrane is in a 
one-phase outside the coexistence curve. 

FIG. 4. Schematic representation of two adhering membranes undergoing a lateral phase separation in the case of the 
bolaform-sticker adhesion. Coexistence between tight (T) and loose (L) membrane domains is shown. The inter-membrane 
distance is smaller than £0 for the tight phase, whereas it larger than £0 for the loose phase. 

FIG. 5. The phase diagram for the bolaform-sticker adhesion when /i_ = as a function of (a) rescaled average composition 
(p+fy^^^ and temperature t/'f, and (b) rescaled average chemical potential ji+fy^^^ and temperature t/^. The continuous line 
is a first-order phase transition line, whereas the dotted line is the limit of metastability of the asymmetric phase with 7^ 0. 
The metastable region of the asymmetric phase is indicated by a thick line in (b). Open circle (o) indicates a critical point. 
Below the critical temperature, there is a coexistence region between the loose (L) and the tight (T) phases as denoted by 
L-l-T. The phase diagram is symmetric with respect to both 0+ — > —0+ and /tt+ — > 

FIG. 6. The phase diagram for the bolaform-sticker adhesion when /x-i- = as a function of (a) rescaled composition 
difference (f>-/'y^^^ and temperature t/^, and (b) rescaled chemical potential difference /x_/7^''^ and temperature t/f. The 
continuous line is a first-order phase transition line, whereas the dashed line is a second-order one. The loose, tight and middle 
phases are denoted as L, T, and M, respectively. The two tricritical points are indicated by a filled circle (•). Below them there 
are two regions of coexistence of the middle phase (0+ = 0) and the tight phase (0+ ^ 0) denoted as Mi-I-T and T-I-M2. The 
phase diagram is symmetric with respect to both </>_ — * —0- and fj, — > — Because /tt+ = there is a degeneracy between 
the tight (T) and the loose (L) phases on the phase diagram. 

FIG. 7. Schematic drawing of the inter-membrane distance of two adhering membranes with coexisting domains between: 
(a) tight (T) and middle (M) phases, and; (b) loose (L) and middle (M) phases, as applies from Fig. 6 for bolaform-sticker 
adhesion. The inter-membrane distance is exactly £0 for the middle phase, whereas it is smaller (larger) than £0 for the tight 
(loose) phase. 



1/3 

FIG. 8. Inter- membrane distance profile at the domain boundary. The variables are rescaled as = Ue S and 
u = Ke^^*Ue^''^x. The values of the two independent parameters are chosen to be J^e/{Kl^^Ue^'') = 0.1 and Ve/Ue^^ = —1. 

FIG. 9. Schematic drawing of the domain boundary and inter-membrane distance of two adhering membranes. Two 
symmetric overshoots in the interface region between the loose and tight domains are shown, in agreement with the model and 
Fig. 8. They result from the combined effect of membrane curvature and lateral tension. 
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